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q ■ ABSTRACT 

5^ | In this paper we present forecasts of the contamination on different shapes of the 

. primordial non-Gaussianity /„/ parameter - detectable on future Cosmic Microwave 

C3 ' Background (CMB) high-resolution anisotropy maps - produced by unresolved ex- 

tragalactic point sources at frequencies of cosmological interest (45-375 GHz). We 
consider two scenarios: an ideal (noiseless) mission and a possible future space-borne 
satellite, with instrumental characteristics similar to the ones proposed for the Cos- 
mic Origins Explorer (COrE). The local, equilateral, orthogonal and fiat shapes are 
t^J- ■ considered in both temperature (intensity) and polarized emission data. The angular 

l/^ | power spectrum and bispectrum of extragalactic point sources are estimated by state- 

of-the-art models of source number counts. The impact of all the most relevant (far-IR 
and radio selected) source populations on these shapes at COrE frequencies is studied. 
\ The results of this analysis show that unresolved extragalactic point sources should 

• not induce a very relevant non-Gaussian signal in the frequency range 100-200 GHz, 

thus not preventing a correct estimate of the CMB primordial f n i parameter. Polariza- 
tion information allows one to significantly reduce the error-bars in the f n i parameter 
and the bias induced by unresolved sources and, hence, to widen the range of frcquen- 
| cies for /„/ studies. On the contrary, at v < 100 GHz or v > 225 GHz, important 

■ non-Gaussian deviations in CMB anisotropy maps are expected due to unresolved 

extragalactic sources. 

Key words: methods: data analysis - cosmic microwave background - extragalactic 
points sources - radio and far-IR: galaxies 



1 INTRODUCTION 

The high precision achieved during the last years to de- 
tect possible weak non-Gaussian signals with primordial 
origin in the most precise CMB anisotropy data maps 
(such as WMAH3 data, the forthcoming Planclfl data 
and future missions) stresses the importance of the char- 
acterization of all the possible non-primordial contami- 
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nants. The standard models of inflation ( Starobinskil ll97^ : 
| Guthlll98ll:lAlbrecht fc Steinhardd 1 1983 : iLindd 1 19821 . 1 19831 : 
iMukhanov et all Il992l ~ predict that the primordial fluctu- 
ations, linearly imprinted in the CMB anisotropics, are 
nearly compatible with a Gaussian distribution. Any de- 
parture from Gaussianity with a primordial origin would 
challenge our understanding of the physics of the early 
Universe. From a practical point of view, the first ap- 
proach to detect non-Gaussian deviations is the study 
of the third order moments (the so-called bispectrum) 
of the anisotropy map. For most of the cases with a 
physical interest, this bispectrum is usually characterised 
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by a single ampli tude, the non-lin e ar coupling param- 
eter f n i (see e.g. Verde et al. 200o]; Komatsu fc Spergell 



120011 ; IBartolo et aJj|2004 iBabich et"aj]|2004lBabichll2005 

For recent reviews of the status of t he primordial 
non-G a ussianity character iz ation s ee e.g. Bartolo et al.l 
1120101); lYadav fc Wandehl d201Ch: iLiguori et al.l (|20ld ); 
IMartfnez-Gonzalez fc Planck Collaboration! (|2013 ). 

The best estimates of the f n i parameter can currently 
be obtained using the publicly available WMAP data. Three 
different shapes with relevant interest for the physics of in- 
flation have been constrained using these data: the local, 
equilateral and orthogonal /„; shapes. Information about 
the inflationary scenarios that can produce these and other 
non-Gaus sian shapes ca n be f o und for example in the re - 
views by IBartolo et ail (|2004h : lYadav fc Wanderd (|201Cj ). 
The instrumental noise contamination and the beam size 
of the WMAP 61 and 94 GHz frequency channels allow 
a maximum resolution of imax ~ 1000. The best esti- 
mates of these shapes provided by the WMAP collaboration 
l|Komatsu et al.ll201ll ) are f^T = 32 ± 21, /** = 26 ± 140 



and f° 



-202 ± 104 (68% CL) for the local equilateral 



and the ort hogonal shapes respectively (see e.g. lCurto et al.l 
l201ll . |2012| . for independent analyses with similar results) . 

Unresolved point sources introduce non-Gaussian 
fluctuations in CMB anisotropics at a detectable level. 
In the past, their "contamination" to the CMB an- 
gular bispectr um has been s t udied in detail by many 



authors (e.g. Refrcgicr et al. 2000; Komatsu & Spcrgel 



200 ll: 



Gorizalcz-Nuevo et al 



lArgiieso et al.l 120031: 
2111)5: lArgiieso et all 120061 : iBabich fc Pierpaolil 



200S 



Serra fc Cooravl I2008T ). As an example, the contamination 



of the f n i paramete r due to point s o urces in WMAP data 
has been studied bv lKomatsu et ail (|2009h considering the 
simple constant-flux model to simulate the radio sources. 
They found that the bias introduced in the local and the 
equilateral shapes is A/^ ; ° c ' ~ 2, A/^ 9 ' ~ 2 2. Taking into 
accou nt the source n u mber counts model by Ide Zotti et al.l 
l|2005h , ICurto et ail (|201lh estimated that the bias in 



fm is A/„' 



(ioc) 



2.5, 



37 and A/£ 



(ort) 



25. 



More recen tly, by explo i ting publicly available all-sky 
simulations, [L et alj (|2012T ) have tried to characterize 
the configuration dependence and the frequency behavior 
of the bispectra originated by "radio" and "far-infrared" 
extragalactic sources. Their results show that the bispec- 
trum due to (extragalactic) far-infrared sources starts to 
dominate that of radio sources on large angular scales at 
150 GHz and dominates on the whole multipole range at 
350 GHz. 

In this paper we characterize the contamination that 
unresolved extragalactic point sources introduce in the CMB 
bispectrum by using the best available predictions on num- 
ber counts of "radio" sources, i.e. Active Galactic Nuclei 
(AGN) with emission spectra dominated by synchrotron ra- 
diation (in this frequency range), and of "far-infrared" (far- 
IR) sources, i.e. low- and high-redshift galaxies in which the 
microwave emission is dominated by warm/cold dust in ther- 
mal equilibrium with the radiation field. We also discuss the 
effect of unresolved polarized extragalactic sources in "fu- 
ture" CMB anisotropy experiments. In particular, we make 
specific predictions for different frequency channels foreseen 
for the future Cosmic Origin Explorer (COrE) satellite mis- 



sion (|The COrE Collaboration et alj|20lj ). which has been 
recently submitted to the European Space Agency (ESA) 
Horizon 2015-2025 Programme^ 

Our current predictions of number counts of extragalac- 
tic sojHxes_are based on two different models: the first one, 
by iTucci et al.l (|201ll ). is able to give the most precise pre- 
dictions on high-frequency (v ^ 100 GHz) number counts 
of radio sources using physically motivated recipes to de- 
scribe the complex spectral behavior of blazars, that dom- 
inate the mm-wave counts a t bright flux densities; the sec- 
ond one, bv lLapi et al.l (|201ll ). using a simple physical model 
and a single reference Spectral Energy Distribution (SED), 
has proved able to reproduce very well all the most recent 
measures (by Herschel, BLAST, etc.) of number counts of 
extragalactic sources at sub-mm wavelengths. It is impor- 
tant to mention that, to our knowledge, this is the first time 
that predictions on the impact of extragalactic radio and far- 
IR/sub-mm sources on the CMB angular power spectrum, 
bispectrum and other statistics are performed with the level 
of detail discussed here. 

In Section [2] we describe the extragalactic point source 
population in this frequency range. The models that char- 
acterise their power spectrum are presented in Section [3] 
Section [4] reviews the basic formalism for the angular bispec- 
trum and the four bispectrum shapes. Section [S] introduces 
the point source bispectra and the estimator to measure the 
A/„; bias. The forecasts for the COrE mission are given in 
Section [S] and, finally, in Section[7]the conclusions are listed. 



2 EXTRAGALACTIC POINT SOURCES 

In this section we describe the expected contribution of ex- 
tragalactic point sources to temperature and polarization 
fluctuations from centimetre to sub-millimetre wavelengths. 
In this range of frequencies, extragalactic sources are typi- 
cally divided in two populations, according to the physical 
radiative mechanism involved: far-infrared (far-IR) sources 
dominated by thermal dust emission and radio sources dom- 
inated by synchrotron emission. In the following analysis we 
use number counts predi c tions obtained by two very recent 
models: i.e., ITucci et alj |201ll . Til, her eafter), as f or ra- 
dio sources at cm to mm wavelengths and lLapi et aL I (|201ll . 
Lll, hereafter), as for far-IR/submm sources. In both cases, 
these models provide predictions on source number counts 
and related statistics which have shown a very good agree- 
ment with observational data. 



3 The COrE mission has not been selected by ESA in the 2011 
first call for proposals. However, a COrE-like space mission with 
technical characteristics similar to the present ones could be hope- 
fully funded in the near future, given that its very high scientific 
output has been already acknowledged by the ESA Space Science 
Advisory Committee (SSAC). In fact, a COrE-like mission will be 
of fundamental scientific importance for detecting the imprint left 
by primordial gravitational waves on polarized CMB anisotropics, 
since the amplitude of primordial gravitational waves is directly 
proportional to the energy scale at which the inflation occurred. 
We don't know the exact instrumental characteristics of this fu- 
ture space mission. However, we do useful predictions by using 
the current instrumental technical det ails quoted in the proposal 
jThe COrE Collaboration et alJliuTlT) . 
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2.1 Radiosources 



2.2 Far-IR Sources 



Recent ev olutionary mod e ls for Extr agalactic Radio Source s 
(ERS) by Ide Zotti etaH (120051 ) and llvlassardi et~al] (|2010l ) 
are capable of providing a good fit to number counts 
from ~ 100 MHz and up to ~ 100 GHz. They adopt 
a simple power-law, i.e. S(u) oc v a , with an almost 
flat spectral index (a ~ —0.1) for extrapolating spec- 
tra of blazar source^] at high frequencies. These mod- 
els, however, over-predict the number counts of extra- 
galactic synchrotron sources detected by the Atacama 
Cosm ology Telescope (ACT) at 148 GHz (|Marriage et all 
l201ll ) and by Planck i n all the High Frequency In- 
strument (HFI) channels JPlanck Collaboration XIIll l 201ll ; 
iPlanck Collaboration Vl3l2013h . The main reason for this 
disagreement is the spectral steepening observed in Planck 
Early Release Compact So urce Catalogue (ERCSC) sources 
at above ~ 70 GHz JPlanck Collaboration VTA l201ll ; 



IPlanck Collaboration xyllgonl) and already suggested by 
other data sets ( Gonzalcz- Nuevo et al 1 l2008l : ISadler et al.l 
I2008D . 

The more recent models by Til use physically grounded 
recipes to describe the complex spectral behaviour of 
blazars, which dominate the mm-wavelengths number 
counts at bright flux densities. The main novelty of these 
models is the statistical prediction of a break frequency, 
um, in spectra of blazars in agreement with classical physi- 
cal mode ls for the synchrotron emi ssion produced in blazar 
jets (e.g.. lBlandford fc Konigll 19791 ) The most successful of 
these models is referred to as C2E:4f|. This model is able to 
give a very good fit to all the data of bright extragalactic 
radio sources available so far: number counts up to ~500- 
600 GHz and on spectral in dex distributions up to, at least , 
200-300 GHz (see Til and IPlanck Collaboration VIll 120131 ; 
lLopez-Caniego et al.ll2012T ). 

Polarization in ERS is typically observed to be a few 
percent of total intensity at cm or mm wavelengths (e.g ., 
iMurphv et aUl2010l ; iBattve et afll201ll ; ISaiina et al1l2011l ). 
with only very few objects showing a fractional polariza- 
tion, n = P/S, as high as ~ 10% . We refer to the recent 
paper by iTucci fe Toffolattil (|2012T ) in order to character- 
ize polarization properties of ERS and their contribution to 
power spectra. In that work authors show that the distri- 
bution of the fractional polarization for blazars is well de- 
scribed by a log-normal distribution with values of the me- 
dian fractional polarization, Timed, and of (II 2 ) 1 ^ 2 depend- 
ing on the frequency. At frequencies v > 40 GHz we take 
n med = 0.04, 0.03, 0.036 and <n 2 ) 1/2 = 0.06, 0.038, 0.045 
for steep-spectrum sources, FSRQs and BLLacs, respec- 
tively. These values are i ntermediate betwe e n the two "ex- 
treme" cases discussed in ITucci fc Toffolattil (I2OI2I) and re ly 
on the most recent data at 86 GHz of Agudo et alj (|2010T ). 



4 Blazars are jet-dominated extragalactic objects, observed 
within a small angle of the jet axis and characterized by a highly 
variable non-thermal synchrotron emission at GHz frequencies, in 
which the beamed compo nent dominates the observed emission 
iAngel fc Stockmarj|l980l'l 

5 It assumes different distributions of break frequencies for 
BL Lac objects and Flat Spectrum Radio Quasars (FSRQs) — with 
the relevant synchrotron emission coming from more compact re- 
gions in jets of the former objects. See Til for more details. 



Recent data acquired by The Herschel Astrophysical Ter - 
aHertz Large Area Survey (H-ATLAS lEales et afl |2010| ). 
the largest a rea survey carried o ut by the Herschel Space 
Observatory (|Pilbratt et aljlioid ). have demonstrated that 
far-IR sources comprise both a low-z galaxy population, 
iden tified through m atchin g to the Sloan Dig ital Sky Sur- 
vey (|York et al.ll2000l) data (|Smith et alJlioill) . and a high-z 
populatio n (median redshift ~2 ) identified through their far- 
IR colors |Amblard et al.l 2010). Low-z galaxies are generally 
nor mal/star-forming late-type galaxies wi th moderate opac- 
ity |Dunne et al.ll201ll ; ISmith et al.ll2012T ). Through analy- 
ses of clustering these two populations are found to be very 
different; the low-z populat ion (z < 0.3) does cluster like 
star-forming blue galaxies (|Maddox et all |2010| ; iGuo et all 
l201ll ; lvan Kampen et alj|2012l ). whereas the high-z popula- 
tion clusters much more strongly, sug gesting that the high -z 
sources reside in more massive halos l|Maddox et al.ll2010D . 

In the case of late-type infrared galaxies, we adopt the 
same source number counts as in Lll, previously derived by 
iNegrello et al.l l|2007h . As explained in more detail in that 
paper, the contributions of the this type of galaxies to the 
number counts were estimated following a phenomenolog- 
ical approach, which consists of simple analytic recipes to 
evolve backwards their local LFs with cosmic time, as well 
as of appropriate templates for their SEDs, to compute K- 
corrections. 

Regarding star-forming proto-spheroidal galaxies, the 
100 and 250 /im luminosity functions (LF) at differ- 
ent redshifts are quite well reproduced by the physical 
model of Early Type Galaxy (ETG ) formation and evolu- 
tion by iGranato et al.l l|200ll . I2004T ), further elaborated by 
lLapi et al. I l|2006t l. and recently revised by Lll without any 
adjustment of the parameters. As discussed in these pa- 
pers, the model is built in the framework of the standard 
hierarchical clustering scenario. This model was the basis 
for the successful predictions of the submil l imetr e counts 
of strongly lensed galaxies by INegrello et all (|2007l ). It also 
accurately reproduced the epoch-dependent galaxy LFs in 
different spectral bands, as well as a variety of relationships 
among photometric, dynamical, and chem ical properties, as 
shown in previous papers (see Table 2 of lLapi et aill2006l ) 
and addit ional result s , espe ci ally on the galaxy chemical evo- 
lution, in lMao et al] l(2007l) ; lLapi et all l|2008l ). 

Number counts from the Lll model are observed to 
be slightly high at mm wavelengthfl The model, in fact, 
tends to overestimate the Poisson contribution to the power 
spectrum of intensity fl uctuations measured at frequencie s 
«: 35 3 GHz by Planc k (IPlanck Collaboration XVIlJ l201lh . 
SPT (|Hall et al.ll2010l ) and ACT (|Dunklev et al.ll201ll ). Scal- 
ing down counts by constant factors is however enough to 
comp ly with measurements at mm wavelengths (|Xia et al.l 
12012? ). In this work we multiply number counts by the fol- 
lowing factors: 0.81 and 0.71 at 375 and 225 GHz, 0.55 for 
1^ 165 GHz. 

Finally, as predicted by iBlainl (1 19961 ): IPerrotta et al.l 



6 This effect is due to the fact that the SED of high redshift 
objects (that contribute to the mm counts around 10-30 mjy) are 
slightly different to the template adopted in Lll as a universal 
SED for protosphcroids. 
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Figure 1. Predicted integral number counts of radio sources 
(solid lines; Til model) and of far-IR sources (total counts: 
dashed lines; proto-spheroidal galaxies only: dotted lines; Lll 
model) at the frequencies indicated inside the panel (each fre- 
quency and the corresponding curves are displayed by a different 
colour). The two black dotted lines, plotted as a reference, repre- 
sent two power laws, with slopes 1 and 3, respectively. 



||2002|. l2003t ) ; INegrello et alj l|2007l ); iPaciga et al] |2009fl ; 
iLima et al.l (|201dl ) among others, the combination of the 
very steep number counts in the far-IR/sub-mm and of 
high redshifts - hence large lensing optical depths - max- 
imizes the fraction of strongly lensed sources. In the case 
of the Lll model, the contribution from lensed proto- 
spheroidal gal axies were estimated using the amplificatio n 
distribution of iPerrotta et al. | (|2003T ): INegrello et all (2007); 
lLapi et al.l (120121). Such p redictions were recently confir med 
by INegrello et al.l (|2010h ; iGonzalez-Nuevo et all (|2012h by 
using the Science Demonstration Phase data of the H- 
ATLAS survey. For the current work we have considered that 
the strong gravitational lensing effect is a random process, 
and therefore, the lensed proto-spheroidal galaxies show a 
Poisson spatial distribution. Due to their lack of correlation 
with the un-lensed proto-spheroidal galaxy sky distribution, 
we have treated them as a different population. 

Regarding polarization, all the - very few - data so 
far collected on the polarization of far-IR/sub-mm sources 
(dusty early- or late-type galaxies and Galactic dust clouds) 
do indicate that the total polarized emission of these sources 
has to be very low with median values at around 1-2% level. 
Moreover, in many observations, no pol arized emission (or 
compatible with ze ro) has been detected |Dotson et al.ll201ol ; 
IVaillancourt]|201ll ). Due to the above, we have assumed - 
for the population of far-IR selected sources as a whole - a 
(conservative) average polarization level of 1%. 



In the case of ERS, we can notice that the slope of integral 
counts (N(> S) oc S~^) is not very different from [3 ~ 1 at 
S < 1 Jy. This means that the contamination of unresolved 
ERS to the CMB power spectra is almost proportional to the 
flux density limit, S c , above which sources can be detected 
and removed or masked. Moreover, because the spectra of 
blazars are typically flat or moderately steep, — 1 < a < 0, 
integral counts change slowly with the frequency, decreasing 
by only a factor 2-4 from 44 to 350 GHz, depending on the 
flux density considered. 

On the other hand, far-IR source counts increase very 
steeply when moving to low flux densities, due to the con- 
tribution of proto-spheroidal galaxies (see dotted lines in 
Fig-HJ. The general slope of far-IR integral counts, /3, is be- 
tween 2 and 3 depending on the flux density range. As a 
consequence, the amplitude of power spectra for unresolved 
far-IR sources is almost insensitive to the value of S c , at least 
for flux density limits typically obtained in large-area CMB 
experiments. Finally, because of the steep frequency spec- 
trum of dust emission, N(> S) is also strongly dependent 
on frequency, making the contribution of far-IR sources al- 
ready completely negligible at frequencies v < 100 GHz (see 
Fig.®. 

In terms of power spectra, it is well known that point 
sources without spatial correlations have: 



Sc 



(I ■> 



(1) 



where k~ x — dBo/dTo is the conv ersion factor from flux 
densi ty to temperature units (see e.g. lTegmark fc Efstathioul 
1996), S c is the flux density limit for the unresolved point 
sources and dn/dS is number counts that characterise 
the population of sources. However, powerful high-z far- 
IR sources have been shown to be highl y correlated by , 
e.g., recent meas ures from He rschel jA mblard et al. 2011; 
Viero et al.ll2012T ) and Planck (jPlanck Collaboration XVIlJ 
201ll ). For the other extragalactic source populations the 



clust ering is negligble (jToffolatti et al.ll2005l : INegrello et al.l 
l2007t ). Therefore, we shall consider the contribution of clus- 
tered sources to angular power spectra (and bispectra) 
only for dusty proto-spheroidal galaxies. This is estimated 
from realistic simulations of t wo-dimensional distribution s 
of clustered point sources (see lGonzalez-Nuevo et al1l2005l ). 
This software was also im plemented in the Planck Sky Model 
jDelabrouille et alj|20l3) and its predictions on the power 
spectrum of cosmic infrared background anisotropies have 
been shown to be in good agreement with the lat est Planck 
measurements l|Planck Collaboration XVIIlll201ll ). 

The total power spectra from point sources are then 
given by the sum of the following components: 



+ c: 



(2) 



3 CONTAMINATION OF POINT SOURCES 
TO CMB POWER SPECTRA 

In Fig. [T] we directly compare expected integrated number 
counts of radio and far-IR sources from the Til and Lll 
models at frequencies ranging from 45 to 375 GHz - i.e. the 
frequencies interesting for CMB studies. The figure well em- 
phasizes the different behavior of the two populations in 
terms of flux density distribution and frequency spectrum. 



where C e and C e are the constant term from 

shot-noise power spectra of radio and far-IR (which includes 
late-type, lensed and proto-spheroidal galaxies) sources, 
and (jsp h ' clust i s the power spectrum from the clustered con- 
tribution of proto-spheroidal galaxiefl 



7 We remind the reader that the total power spectrum of a pop- 
ulation of clustered sources is C\ ot = Cf hlte + Cf u3t . 
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Figure 2. Point source power spectra for temperature (top panels) and polarization (bottom panels) at 75, 105, 165 and 225 GHz. Lines 
in colour are for radio sources (solid lines) and for far-IR sources (dashed lines for the total contribution; dotted lines for the Poisson 
contribution). We consider two flux density limits for radio sources: S c =1 and 0.3 Jy; for far-IR sources we plot only the case with 
S c =1 (there are not significant differences taking S c = 0.3). Black lines are for CMB spectra from the fiducial WMAP 7-year model: in 
the bottom panels dotted lines are for the E-mode, dashed lines for lensing-induced B-mode and solid lines for primordial B-mode with 
r = 0.1, 0.01, 10" 3 . 



Regarding polarization, point sources contribute, on av- 
erage, equally to power spectra of E- and B-modes or equiv- 
alently to Q and U power spectra in the flat-sky limit. Po- 
larization power spectra of point s ources are propor tional to 
the temperature power spectrum dTucci et al.ll2004h : for the 
Q power spectrum, e.g., we have 

C? = kl (S 2 p 2 cos 2 (20)) = ^ (p 2 ) Cf . (3) 

The above relation assumes that the flux density S, the 
polarization fraction p and the polarization angle in the 
chosen reference system <\> are independent variables (see 
iTucci fc Toffolattill2012l , for a detailed discussion on this as- 
sumption). On the other hand, all the cross-power spectra 
are equal to zero: 

Cj Q oc (S 2 p cos(2(t>)) = 
Cf u oc (S 2 p sin(2</>)) = 

Cf v oc (S 2 p 2 cos(20) sin(20)> = . (4) 

Fig. [2] shows model predictions for temperature and polar- 
ization power spectra at 75, 105, 165 and 225 GHz. We also 
consider two possible flux density limits, i.e. S c = 1 and 
0.3 Jy. 



As previously discussed, far-IR spectra are almost in- 
sensitive to the choice of S c , especially for the clustered 
contribution. Comparing the two source populations, we see 
that radio sources are dominant at v < 100 GHz, while far- 
IR sources start to dominate at v > 200 GHz both in tem- 
perature and polarization. At intermediate frequencies, i.e. 
100 < v < 200 GHz, ERS are the main contaminant at high 
multipoles (£ > 1000), but not for I < 1000, where the far- 
IR clustered sources contribution becomes higher than the 
ERS power spectrum in temperature and comparable in po- 
larization. 

Because proto-spheroidal galaxies are the only clus- 
tered population in the Lll model, it is worth discussing 
more in detail their relevance in terms of power spec- 
trum. In Fig. [3] we plot the ratio of the total spectrum 
of proto-spheroidal galaxies to their shot-noise spectrum, 
Qsph ^£,s P h,white ^ ^ shot-noise spectrum of all far- 
IR sources, C^ ph /C l e r ' whlte , (the former ratio enters later 
in Eq.[TS] and gives us the increase of the bispectrum due 
to clustering). As the frequency decreases, we observe two 
effects: o) the contribution of proto-spheroidal galaxies to 
the far-IR power spectrum becomes less and less impor- 
tant; b) moreover, the relative contribution of clustering 
to proto-spheroidal power spectrum reduces. The reason of 
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Figure 3. The ratio of the total spectrum, C^ ph , to the shot- 
noise spectrum, c'P h ' whlte f Q f p ro to-spheroidal galaxies (dashed 
lines), and to the shot-noise spectrum from far-IR sources, 
gir.white. ^gQijjj i mes ^ a t the frequencies indicated in the panel. 



these trends is related to the shift of the redshift distribu- 
tion of submillimetre bright galaxies to higher redshift with 
increasing the wavelength, as effect of the strongly nega- 
tive A'-correction (Lll). At high z, however, massive proto- 
spheroidal galaxies become increasingly rare and the slope 
of number counts at low frequencies is not as steep as that 
at sub-mm wavelengths. As a consequence, the shot-noise 
term of the power spectrum, dominated by bright sources, 
has more weight with respect to the clustered one (domi- 
nated by faint sources). 

Finally, we discuss the contamination of extragalactic 
sources to CMB fluctuations: in the frequency range consid- 
ered, it is relevant only at very high I (> 2000) for tempera- 
ture and is essentially negligible for the E-mode polarization 
(see Fig. [2]). This is not the case for the B-mode polariza- 
tion arisen from primordial gravitational waves if r < 0.01. 
A subtraction of points sources at levels even lower than ones 
shown in Fig.[2] is not an easy task. For radio sources, this 
could be in principle done a) by removing/masking sources 
down to flux density limits of S < 0.3 Jy or b) by a direct 
measurement of the power spectrum of polarization fluctua- 
tions at small angular scales. In both cases observations with 
arcmin resolution are required. Far-IR sources are still more 
challenging to deal with, and the previous methods are not 
valid to clean their contamination due to the high correlated 
signal. Instead, approaches similar to ones used for diffuse 
Galactic foregrounds should be more suitable for this goal, 
although they are model dependent. 



4 THE ANGULAR BISPECTRUM 

The CMB angular bispectrum is defined in terms of cubic 
combinations of the spherical harmonic coefficients of the 
CMB anisotropies of the T, E and B modes: 



k 



(5) 



where the indices {i,j,k} cover the T, E and B modes. 
The angular- averaged bispectrum and the reduced angular- 



averaged bispectrum, B l J k , „ and b l j e . are defined as: 













12 h 


m 1 m 2 m 3 \. 


mi 


7712 7773 


ijk ( $-1 


4 




I o 








(6) 



The reduced bispectru m can be written as (see, e.g. 



lYadav fc Wandel 



spectru 
t 2010) 



(47t) 3 H)* 1+ ^ 3 



mi 



«2 «3 
T772 7773 



J (2tt) 3 (27) 3 (2tt) 3 



z k 2 d 3 k 3 



(ki)Y e * 2m2 (k 2 )Y e * 3m3 (kz) 



<7l 1 (fci)5i(fe2)5 < fc 3 (fc3){$(k 1 )$(k 2 )<I>(k3)>, 



where gl(k) is the transfer function for temperature or polar- 
izatior|f| and (<I>(ki)<l>(k2)$(k3)) is the 3-point correlation 
function in Fourier space of the Bardeen's curvature pertur- 
bations ( I > (x). The 3-point correlation function can be given 
in terms of the so-called shape function F(ki, k 2 , fe): 

($(ki)$(k 2 )$(k 3 )) = (2 7 r) 3 5 3 (k 1 + k 2 + ka)F(fci, k 2 , k s ) 

(8) 

The level of primordial non-Gaussianity is parametrised 
by t he non-linear couplin g para m eter /„; dVerde et all 
20001: iKomatsu fc Spergell 1200 ll : iBartolo et alj 120041 : 



Yadav fc WandeltJ |2010). This parameter measures de- 



partures from zero of the bispectrum. Depending on the 
physical mechanisms of the different inflationary models the 
shape function can take different forms. The shapes that we 
consider in this work are the local, equilateral, orthogonal 
and flat shapes. 

• Local shape. This type of Non-Gaus sianity is gen- 
erate d in multi - field inflationary models (Komats u et all 
120051 : iKomatsul I2010T ). the curvaton dLvth et all I2003TI 



the inhomogeneous reheating scenario (IDvali et al 
IBartolo et ail |2004 ), hybrid inflation 

lYadav fc WandeltfeOlJ ) 



Its shape is given by (see e.g. ICreminelli et al 



(iLinl il2009h . 



2004; 



etc . 
120061 : 



F(k u k 2 ,k 3 ) = 2A 



'fnl [ 



3-(n s 



-1)7. 3 - 

/t 2 



— + (2 perm)] , 



where A is the amplitude ot tne power spectrun 
and n s is the spectral index (|Komatsu et alj 1201 ll ). 
Most of the signal of this type of bispectrum is 

located i n squeezed configura t ions such as fc i ^ 

k 2 » fca llCreminelli et all 120061: lYadav fc Wandeltl I2OI0I : 



the amplitude of the 



(9) 
ectrum 



iMartmez-Gonzalez fc Planck Collaboration! 20121 ) 

• Equilateral shape. This type of Non 
Gaussian ity is generated in the Dirac-Born-Inf eld in- 



flation (ISilverstein fc Tond 120041 : IBartolo et al] 12004 



Langl ois et al.ll2008l ). ghost inflation (|Arkani-Hamed et al.l 
2004h . single-field inflationary models in Einstein gravity 



8 This function ca n be computed numerically using the CAMB 
llLewis et al.ll200oh software |http: / / camb.info| 
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(IChen et al, 2007) et c. Its shape is given by (see e.g. 
00d; r 



ICreminelli et alj|200d lYadav fc Wandeltj|2010h 

1 



F(ki,k2,k 3 ) = 6A'f n 



3_(„ s _l) 3_(„ s _l) 



Yy + (2 perm) 



+ 



(fclfc 2 fc 3 ) 2 ( 4 -™=)/ 3 
1 



{fc< 4 



n.)/3 fc 2(4-»i,)/8 i .(4 



— + (5 perm)}]. 

(10) 



Most of the signal of this type of bispectrum is located 

in equila teral configurations su c h as fci ~ fc 2 ^ fcj 

(see e.g. [Creminelli et all l200d: lYadav fc Wandeltl l20ld ; 



iMartmez-Gonzalez fc Planck Collaboration! 20121 ). 

• Orthogonal shape. This type of N on-Gaussianity 
is ge ne rated in gener al sing le-field models l|Cheung et al.l 
120081 ; ISenatore et al.l l20ld). I ts shape is given by 



_ .. JJ . i ts snape is given Dy 
(see e.g. ISenatore et all l20ld ; lYadav fc Wandeltl l20ld ; 
iKomatsu et al.ll201ll ): 



F(k!,k2,k 3 ) =6A 2 / ni [- 



3-(n s -l)j,3-(n 3 -l) 



Yy + (2 perm) 



+ 



(fcifc 2 fc 3 ) 2 ( 4 -™=)/ 3 
( 3 



, ^-„,)/3, 2(4-ra a )/3, (4-n s ) 
r,j ft, 2 ^3 



+ 



(5 perm) J-l . 

(11) 



This bispectrum peaks in both equilateral fci ~ fc 2 — 
&3 and flat-triangle configurations such as fci = fc 2 + 
fc 3 (see the shape figures e.g. in ISenatore et ah! l20ld ; 



IMartmez-Gonzalez fc Planck Collaboration! 20121 ). 

• Flat shape. This type of Non-Gaussianity is gener- 
ated for example in models with modifications in the i nitial 
state of the inflaton field (see e.g. iMeerburg et alj|2009l . and 
references therein) , single- field models with high derivative 
iterations jBartolo et al.ll2010bh among others. Its shape is 
given by (|Meerburg et al.ll2009l ): 



F{k 1 ,k 2 ,k 3 ) = 6A 2 f nl [ 



3-(n a -l)i,3-(n s -l) 



YY + (2 perm) 



+ 



(fcifc 2 fc 3 ) 2(4 -" s) 



+ 



(5 perm) | 

(12) 



This bispectrum peaks in flat- triangle configurations fci 
fc2 + fc3 (see the shape figure in IMeerburg et alj|2008h . 



Using the shape forms defined in EQS.I91 fT0lllll . [T2l and the 
CMB bispectrum in Eq. [7| we have: 

fe* 3 ) ^ = 2 J x 2 dx[a^ (x)pi {x)pt 3 (x) + (2 perm)] , 

(13) 

/DO 
dxx 2 y - ot\ (x)P 3 l2 (x)Pt 3 (x) + (2 perm) 

( x hi 2 (x)Se 3 (x) + (5 perm) - 26l t (x)S J i2 (x)Se 3 (as)] , 

(14) 



KUr* = 18 \ n dxx 2 [ - ^[x)pi{x)pl{x) + (2 perm) 

(15)" 



+Pi 1 (x) 7 i 2 (x)Sl(x) + (5 perm) - |4 (^(x)^ (*)] 



and 



(%" 2 ^ 3 ) !J = 6 / ctea; al 1 (a;)^ 2 (a;)ft 3 (a;) + (2perm) 



h flat _ g / , . J I , 

(^h^M^M + (5 perm) + 3<5^ (x)5j 2 (aO$/ 3 (a;)] 

(16) 



where a\(x), P\(x), 7l(x), 5|(:r) are filter functions defined in 
terms of the po wer spectrum P(k) and th e transfer functions 
q}(k) (see e.g. iFergusson et al.l [2010al lbl lYadav fc Wandeltl 
l20ld ; ICurto et al. 2011 , for their actual definitions) . 



5 BIAS AND UNCERTAINTIES ON 
PRIMORDIAL NON-GAUSSIANITY 

5.1 Bispectrum of point sources 

The bispectrum of a point source population without spatial 
correlations, as the case of radio sources, is well characterised 
and is constant with I: 



i white 



„3 dn 

S ds dS - 



(17) 



On the other hand, far-IR sources have spati al correlations. 
Accor ding to the widely used prescription bvlArgiieso et ail 
e.g. in lGonzalez-Nuevo et al.l (|2005f ); 
and extended to the full sky by 



Argiicso et al. 


(hood) 


Lacasa et al.l ( 


2012), t 



tered sources population can be written as: 



iwhite 
Ot 1 l 2 l 3 



r'tot rnot r'tot 
°<?i U t 2 U l 3 

/~i white /~i white white 



(18) 



The total bispectrum of point sources is then given by the 
sum of the shot-noise contributions of ERS (^/ 2 '™ 3 '" te ) and 
of far-IR late- type and lensed galaxies (b l ^ 2 e e ' wh ' te ) plus the 
"clustered" bispectrum of proto-spheroidal galaxies (6|^ 2 ^ ) 
as given by Eo. 1181 



h ps 

°*1*2*3 



_rrad, white . ilt-\-le,white . isph 

-°e 1 t 2 e 3 + eit 2 i 3 + °eit 2 £ 3 



(19) 



When polarization is considered, we have 10 different bis- 
pectra arisen from all the possible combinations of temper- 
ature and polarization E— and B-modes. However, on the 
base of the discussion in Section[3] it is easy to find out that 
the cross-bispectrum bJ^ 2 B i3 and all bispectra with an odd 
number of polarization terms are zero, i.e. bf E 2 E t3 = bj^ 2 £ = 
... — 0. For the remaining two cases, the bispectrum is pro- 
portional to bj^ &3 and it is given as: 



iTEE _ iTBB _ 1 / 2\,TTT 
°l^l 2 l 3 — °*1*2«3 _ 2 °^2«3 



(20) 



The total power spectrum and the bispectrum of the point 
sources are the two ingredients needed to provide an esti- 
mate of their impact on the primordial non-Gaussianity /„; 
parameter. The shot-noise terms of TTT and TEE bispec- 
tra from radio and far-IR sources are shown in Fig. [4] The 
frequency at which far-IR sources equal or overcome the 
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Figure 4. The amplitude of the shot— noise bispectrum of radio sources (b rad ) and of the far-IR source (for late-type plus lensed galaxies, 
b lt + le , and for proto-spheroidal galaxies, f, s P' l > ,i, '" t ' : ) as a function of frequency and for the two considered flux density limits, i.e. Sc = 0.3 
and 1 Jy. The left panel corresponds to the TTT bispectrum and the right panel corresponds to the TEE bispectrum. 



local 




equilatera 



100 



000 



max multipole I 
orthogonal 





max multipole 



Figure 5. Ratio of the far-IR source bispectrum b lr to the radio source bispectrum b rad for S c = 0.3 Jy and the 4 considered shapes: 
local configuration (£\ = £2 = £ and £z = £ m in = 2), equilateral configuration (£\ = £2 = £3 = £), isosceles orthogonal configuration 
(£± = £2 = £ and £3 = \/2£) and isosceles flat configuration (£\ = £2 = £ and £3 = It). The solid lines correspond to T, and the long 
dashed lines correspond to T and E. 



ERS bispectrum is shifted to higher frequencies with respect 
to the power spectrum, i.e. it is around 300 (350) GHz for the 
TTT (TEE) bispectrum. This is in agreement with the fact 
that sources with flux density close to S c give a larger con- 
tribution to the bispectrum than to the power spectrum. A 
significant dependence on the flux density limit is observed 
in the bispectrum of ERS and of late-type plus lensed galax- 
ies, but not of proto-spheroidal galaxies. 

Finally, we compare the total bispectrum of far-IR 
sources with the ERS one as a function of multipoles and 
for the four used configurations (see Fig.0. The ratio is > 1 
only at 225 GHz for TTT, and for TEE, but at low mul- 
tipoles. Otherwise, far-IR sources are always subdominant, 



or even negligible, unless considering multipoles I < 100. As 
widely discussed in Sec.O clustering of far-IR sources be- 
comes less relevant at low frequencies, and at v < 100 GHz 
the far-IR bispectrum is almost constant. Our Fig . [5] can 
be directly compared to Fig. 9 of lLacasa et al.l (|2012h . Var- 
ious differences can be highlighted from the two plots: 1) 
at high multipoles, the bispectrum ratio we find is typically 
higher than one or der of magnitude; 2) far-IR bispectra by 
lLacasa et al.l (|2012h have a steeper slope; 3) n o frequency de- 
pendence is observed in|L et al.1 (|2012l ) results. These 
discrepancies should be attributed to the different way to 
model far-IR sources and their correlation. 
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5.2 Bias and uncertainties on /„; 

Considering weak levels of non-Gaussianity, the bias induced 
in the primordial bispectrum B ?!'™ by the point-source 
bispectrum B?\ , is (see e.g. Komatsu fc Spergell l200ll: 
Bartolo et alj BoM IVadav fe Wandeltl l20ld; lLacasa et all 



20121 ; iMartmez-Gonzalez fc Planck Collaborationll2012l ) 



where 



n e 1 e 2 e 3 - D e 1 £ 2 e 3 



r*tot /~>tot s~itot\ 



(21) 



Of =(0(7 + G£ juv+Gf , (22) 

Cf MB is the CMB power spectrum, contains the in- 
strumental and pixelization beams, C'" mae is the power 
spectrum of the noise, A tl i 2 i 3 = 1 + 25e 1 e 2 Si 2 e 3 + &t x t 2 + 
+ 5e, i 3 . The error-ba r of the f n i parameter is given by 



ijKomatsu fc Spergell lioOlh 



\- D t 1 e 2 t 3 ) 



r*tot c*tot ritot\ 
i^2^3 K—tiMa^ii U l 2 ) 



(23) 



This estimator can be generalised to include the polariza- 
tion: 



A/ni = cr 2 (/ n! )x 

ijkrat ii^l 2 ^:i 3 ^i 



,; ■st(-Sf 1 f 2 < 3 ) > (24) 



where 



{!X 



1 ijk 



J e 1 e 2 e 3 J 



CTlrAB^Zj 3 ( 25 ) 



ijkrat £ l ^2 ^3 =S^maic 

The covariance matrix Cijk, rs t is given by (|Yadav et al.l 
120071 ): 

<^ijk,prq — ^£ 2 <^l 3 + £l <^i 2 O j>3 0t 2 l 3 + Ls £i l^i 3 0l x i 2 

+ c t r 1 C J e ^C^5e 1 e 3 + Cl q 1 C£C^5 ei e 2 5e 2 e 3 5 e3 e 1 

+ c l r i C Z c 'i 3 s iii3 5 t 2 e 1 5e 2 i 3 . (26) 

The indices {i, j, k, r, s, t} cover the intensity T and the po- 
larization E and B modes. The cosmological signal for the 
B mode is negligible. Including the point source B mode sig- 
nal in the estimator does not improve the estimates of cr(/„;) 
and A (/„;). Therefore we just include in the estimator the 
T and E modes for the case with polarization. As in the 
case with intensity only, each of the power or cross-power 
spectrum in Eq. ([26} contain the CMB signal, instrumen- 
tal beam, instrumental noise spectrum and the point source 
spectrum contributions. 



6 RESULTS 

We have studied the impact of unresolved point sources on 
the possible measurement of the f n i parameter for two cases: 
a) an ideal experiment without instrumental beam nor in- 
strumental noise; b) a possible future COrE-like mission, 
sensitive to temperature an d polarization, with instrumen- 
tal parameters adopted from lThe COrE Collaboration et al.l 



Table 1. Instrumental properties for a possible future COrE-like 
mission, based on values given by The COrE Collaboration et al. 

toil) . 



v (GHz) 


45 


75 


105 


165 


225 


375 


FWHM (arcmin) 


23.3 


14.0 


10.0 


6.4 


4.7 


2.8 


T RMS noise* 


5.25 


2.73 


2.68 


2.67 


2.64 


68.6 


Pol. RMS noise* 


9.07 


4.72 


4.63 


4.61 


4.57 


119 



* fiK arcmin 



(2011) and summarised in Table[T] The CMB power spec- 
trum and transfer function for temperature and polarization 
needed to estimate the f n i bias and its uncertainty (Eqs. 
[22 rj31 EH a nd 1231) are numerically computed with CAMB 
I Lewis et al. 2000) using the cosmological parameters that 
best-fits WMAP 7-yr data. In this section we present results 
for the uncertainty on the /„; parameter, a(f n i), for the ex- 
pected bias, A(/ n ;), and for the relative bias, A(/ n ;)/cr(/ n ;) 
(see Fig.[6] Tables[2] and [3] for the ideal case, and Fig.0 Ta- 
ta les[H and [S] for a COrE-like mission). We compare results 
obtained with temperature only and with temperature plus 
polarization. Four different f n i shapes are taken into ac- 
count, i.e. local, equilateral, orthogonal and flat. Fig. [S] and 
[7] are for S c = 0.3 Jy and for frequencies 75, 105, 165 and 
225 GHz. In the Tables we also provide the results, evaluated 
at £ max = 2048, for S c = 1.0 Jy and for 45 and 375 GHz. 

The following results assume full sky. This is an unreal- 
istic hypothesis as Galactic Plane and strong point sources 
have to be masked. How cr(/„;) depends on the mask and on 
the fraction of the sky available, f s k y , is something complex 
and outside the purpose of this work. However, for large ar- 
eas of th e sky (f s kv > 0-7) the following approximation can 
be used i|Komatsu et al .1120021 '): 

£>"(/n()cutsky = O {fni ) fullsky / sj fsky ■ (27) 

The bias on /„; does not depend on fsky, as it can be seen 
from Eq.EU 



6.1 Uncertainty 

In Eqs. l23l and l25l of the previous section we have provided 
the analytical formulas to estimate the uncertainty on the 
fni parameter. As we can see, the point sources contribution 
enters only in the total angular power spectrum (Eq. l22p . in 
the sum Cf MB + Cf . In the frequency range 75-225 GHz, 
C\" is always subdominant with respect to T and E-mode 
CMB spectra up to £ ~ 2000 (see Fig.rjJ. Therefore, as il- 
lustrated by Fig. [6] and [7] the impact of point sources on 
the uncertainty in the determination of the /„; parameter 
is small or negligible, especially when polarization is consid- 
ered. In this case, in fact, a{f n i) is almost indistinguishable 
with the intrinsic f n i uncertainty due to the cosmic vari- 
ance (black curves in the plots, indicated as "no sources"). 
This is true also for the "temperature only" case at 105 and 
165 GHz, i.e. the frequencies where the signal from point 
sources is minimum. At 75 GHz the uncertainty becomes 
slightly higher only at £ > 1500 because of the increasing 
contribution of ERS; while, a larger a(f„i) is obtained at 
225 GHz in the all range of is, clearly due to the clustered 
signal of IR sources. At 353 GHz the point sources contam- 
ination grows dramatically, and the <j(fni) is around 4-6 
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Table 2. The expected uncertainty cr(/ n ,), bias A/ n , and relative bias A/ n ;/(r(/ n ;) at l m ax = 2048 from temperature only and 
temperature plus polarization for an ideal mission. S c = 0.3 Jy is adopted. 





case 


T 


T 


T 


T 


T 


T 


T+E 


T+E 


T+E 


T+E 


T+E 


T+E 




Frcq. 


45 


75 


105 


165 


225 


375 


45 


75 


105 


165 


225 


375 




v(fnl) 


7.4 


4.9 


4.4 


4.2 


5.2 


29.7 


2.8 


2.3 


2.1 


2.1 


2.3 


3.3 


local 




82.1 


10.3 


2.1 


0.4 


0.7 


25.6 


-9.7 


-2.6 


-0.6 


-0.1 


-0.1 


-0.1 




A/„,/<r(/ n ,) 


11.1 


2.1 


0.5 


0.1 


0.1 


0.9 


-3.5 


-1.1 


-0.3 


-0.1 


-0.1 


-0.0 




<r(fnl) 


75.4 


60.4 


56.8 


55.9 


65.1 


232.4 


24.0 


20.3 


19.5 


19.2 


21.1 


35.3 


equilateral 


A/ n , 


4629.3 


599.9 


107.2 


17.2 


49.1 


950.7 


480.2 


54.5 


7.3 


0.8 


4.0 


17.9 




Af nl /a(f nl ) 


61.4 


9.9 


1.9 


0.3 


0.8 


4.1 


20.0 


2.7 


0.4 


0.0 


0.2 


0.5 






40.1 


29.7 


27.4 


26.8 


32.0 


139.0 


11.6 


9.8 


9.4 


9.2 


10.1 


16.5 


orthogonal 


A/ n , 


955.5 


51.7 


2.5 


-0.2 


9.7 


371.4 


308.7 


48.3 


9.4 


1.6 


3.4 


7.2 




A/„,M/„z) 


23.8 


1.7 


0.1 


-0.0 


0.3 


2.7 


26.7 


4.9 


1.0 


0.2 


0.3 


0.4 




""(/nj) 


73.7 


53.3 


49.0 


47.9 


57.6 


284.7 


20.6 


17.4 


16.6 


16.4 


17.9 


29.1 


flat 


A/„, 


598.6 


150.5 


36.0 


6.7 


3.4 


-66.0 


-311.6 


-56.3 


-12.2 


-2.2 


-3.9 


-5.2 




A/ n ,/<r(/ n ,) 


8.1 


2.8 


0.7 


0.1 


0.1 


-0.2 


-15.1 


-3.2 


-0.7 


-0.1 


-0.2 


-0.2 



Table 3. The expected uncertainty <r(/ n ;), bias A/ n , and relative bias A/ n ;/(r(/ n ;) at lmax = 2048 from temperature only and 
temperature plus polarization for an ideal mission. S c = 1 Jy is adopted. 





case 


T 


T 


T 


T 


T 


T 


T+E 


T+E 


T+E 


T+E 


T+E 


T+E 




Freq. 


45 


75 


105 


165 


225 


375 


45 


75 


105 


165 


225 


375 




0"(/ni) 


9.9 


5.9 


4.9 


4.4 


5.3 


29.8 


3.2 


2.5 


2.3 


2.1 


2.3 


3.3 


local 


A/n, 


351.1 


68.5 


17.6 


3.4 


2.3 


25.6 


-5.9 


-12.6 


-4.4 


-1.0 


-0.5 


-0.1 




A/n,/<7(/n,) 


35.6 


11.6 


3.6 


0.8 


0.1 


0.9 


-1.9 


-5.1 


-2.0 


-0.5 


-0.2 


-0.0 




0"(/ni) 


88.7 


66.6 


59.9 


57.0 


65.6 


229.8 


27.7 


21.8 


20.2 


19.5 


21.2 


35.3 


equilateral 


A/n, 


17689.6 


4075.1 


1017.5 


174.2 


135.5 


948.9 


2226.3 


390.7 


90.7 


11.5 


10.8 


20.5 




Af nl /a(f nl ) 


199.4 


61.2 


17.0 


3.1 


2.1 


4.1 


80.4 


17.9 


4.5 


0.6 


0.5 


0.6 




f(/ni) 


50.2 


33.8 


29.3 


27.5 


32.3 


138.6 


13.4 


10.5 


9.7 


9.4 


10.1 


16.5 


orthogonal 


A/n, 


5411.3 


573.7 


80.9 


3.6 


16.4 


377.8 


1007.2 


297.6 


82.1 


15.1 


9.3 


7.1 




A/n,Ar(/„,) 


107.8 


17.0 


2.8 


0.1 


0.5 


2.7 


75.3 


28.4 


8.5 


1.6 


0.9 


0.4 




f(/ni) 


94.6 


61.3 


52.7 


49.2 


58.2 


284.5 


23.8 


18.6 


17.3 


16.6 


18.0 


29.2 


flat 


A/n, 


455.7 


784.2 


263.4 


59.1 


26.6 


-69.0 


-770.4 


-326.8 


-96.4 


-19.6 


-10.8 


-4.1 




A/ n ,/cr(/ n ,) 


4.8 


12.8 


5.0 


1.2 


0.5 


-0.2 


-32.4 


-17.5 


-5.6 


-1.2 


-0.6 


-0.1 



times higher than at 225 GHz. On the other hand, the loss 
of sensitivity between 75 and 45 GHz is not so strong, due to 
the natter spectral behaviour of ERS. These results seem to 
be quite independent of the two flux density limits consid- 
ered, even at 45 and 75 GHz, where a(f n i) is about 10-20% 
larger (see Tabled}©. 

When polarization is taken into account, / n , error-bars 
are reduced approximately by a factor of 2-3 accordingly to 
the shape, in the frequency range 75-225 GHz. At 45 and 
375 GHz, polarization permits to achieve a very large im- 
provement in the / n , error-bars, but only for the ideal case. 
In a COrE-like mission the instrumental noise in polariza- 
tion saturates the f„i error-bars at these frequencies and 
prevents to reduce them as in the ideal case. 

The values of cr(f n i) obtained from a COrE-like mission 
are very close to the results of the ideal case. This means 
that observational performances, in terms of instrumental 
noise and beam, given in Table[T] are almost optimal for f n i 
measurements, at least up to i ~ 2000. Only in the 75-GHz 
channel (and partially in the 105-GHz channel) the effect of 
the beam is observed in Fig.[7] where a(f n i) starts to be con- 
stant at lmax > 1500. In the 165- and 225-GHz channels the 



/„, uncertainty is always decreasing up to £ ma x ~ 2048. Note 
however that this is not the case in a Planck-like mission, 
where the noise contamination in the E mode is more impor- 
tant and the i mprovement from T only to T and polariza- 
tion is 20-30% (M artinez-Gonzalez fc Planck Collaboration! 
I2012D . 



6.2 Bias 

As expected, the bias (or the relative bias, shown in Fig. [6] 
and [7]) rapidly increases with i m ax- Its value is also strongly 
dependent on the frequency. In fact, for example in the case 
of no clustering, the bias is directly proportional to the am- 
plitude of the shot-noise bispectrum of point sources (see 
Eq.IUll- 

If we consider only temperature and S c = 0.3 Jy, the 
bias is practically negligible at 165 and 225 GHz (i.e. where 
b^"' e ^^ te is minimum) for the local, orthogonal and flat 
shapes. The relative bias is less than 1 also for the equilateral 
shape, which is the shape most affected by point sources. 
As for a(f n i), the bias strongly increases at 375 GHz, be- 
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ing typically more than one order of magnitude higher than 
at 225 GHz. At 105 GHz it is still small, with a relative 
value less than 1 (apart from the equilateral shape, with 
A/ni = 107). At lower frequencies, however, A/ n ; 3> 0-(/„;). 

Moreover, we find that A/ n ; for the local, equilateral 
and flat shapes is negative, although negligible, at low mul- 
tipoles (£ ^ 200) and it becomes positive at higher mul- 
tipole^j]. The bias for the orthogonal shape is oscillating 
around zero and therefore, in this case, there are certain 
multipoles where Af n i ~ 0. 

It is interesting to notice that when instrumental noise 
and beam function are introduced, the bias (and the relative 
bias) are equal or even lower than in the ideal case at £ m ax = 
2048. One simple explanation is that the beam and the noise 
prevent the high-^ terms in the sum of Eq.[2U and [2H to 
diverge. 

When the polarization is considered, the bias suffers im- 
portant variations for each different shape in its amplitude 
and its sign. This is dependent on the extra terms in the 
estimator: the E-mode primordial bispectrum, the E-mode 
point source bispectrum and also the inverse of the covari- 
ance matrix in Eq. 1261 The bias of the local and flat shapes 
becomes negative, whereas in the equilateral and orthogonal 
shapes it keeps a positive value. The local and flat shapes 
bias are also reduced in amplitude but the strongest im- 
provement is for the equilateral shape. On the contrary, at 
certain frequencies the amplitude of the bias of the orthogo- 
nal shape increases compared with the case of temperature 
only (see Tables [2] and [3} . Finally, it can be noticed that 
the amplitude of the bias is typically lower in a COrE-like 
mission than in the ideal effect of considering in- 

strumental properties in Eq.l21land l24l 



7 CONCLUSIONS 

In this paper we present a full set of forecasts on the contam- 
ination produced by unresolved extragalactic point sources 
on the primordial non-Gaussianity /„; parameter - hope- 
fully detectable by the analysis of future high-resolution 
CMB (all-sky and low-noise) anisotropy maps - in the fre- 
quency interval 45-375 GHz. We have considered two sce- 
narios: a) an ideal case without instrumental noise nor in- 
strumental beams; b) a future experiment with COrE-like 
instrumental parameters. We have characterized all source 
populations relevant in this fr equency interval by using two 
recently published models: the lTucci et al.l (12011) m odel for 
extragalactic radio sources and the Lapi et alj 1 201l[ ) model 
for far-IR/sub-mm sources. It is important to remind here 
that these two models have proved very successful in the 
explanation of all of the most recent and published results 
on number counts (and on other relevant statistics) of ex- 
tragalactic point sources in the frequency range from ~ 30 
to ~ 1000 GHz. 

As a first step, the angular power spectrum and bis- 
pectrum due to unresolved extragalactic point sources for 
both temperature and polarization have been estimated and 



9 This change in the sign of the "local" bias is not found in 
lLacasa et al.l ||2012T ), where the bias due to IR sources is still neg- 
ative at I = 700 and 2048. This different result may be dependent 
on the actual details of the characterization of the IR sources. 



compared. According to the above quoted models, the con- 
taminating signal produced by far-IR/sub-mm selected ex- 
tragalactic sources starts to dominate over the ERS one at 
v > 200 GHz in the temperature (or total intensity) angular 
power spectrum and bispectrum. In spite of their very low 
(average) level of polarized emission, this source population 
is also relevant in polarization, especially at large scales, 
due to the enhanced signal produced by the strong clus- 
tering of high-z far-IR/sub-m spheroidal galaxies. On the 
contrary, their effect is completely negligible at frequencies 
< 100 GHz. On the other hand, unresolved ERS constitute 
a more relevant contaminant in all the frequency range con- 
sidered here and, thus, their effect on the CMB anisotropy 
signal has to be properly modeled, for a correct evaluation 
of their impact. In any case, their signal is found to be dom- 
inant over the intrinsic CMB one only at frequencies below 
75 GHz, thus out of the CMB cosmological window. 

We have studied the expected impact of undetected 
point sources on the /„; parameter in terms of uncertainty 
and bias, for the local, equilateral, orthogonal and flat 
shapes. For a COrE-like mission, we find out that the best 
frequency for a f„i analysis is ~ 165 GHz. At this frequency 
(and S c ~ 0.3 Jy), the effect of unresolved point sources on 
the accuracy in the determination of f n i is almost negligi- 
ble. The uncertainty on /„; does not significantly change 
whether or not their contribution is included, and estimated 
values of /„; are essentially unbiased (apart from the equilat- 
eral shape for which Af„i/a(f n i) ~ 0.3 at £ ma x = 2048). At 
the frequencies 75 and 105 GHz, results are similar to 165- 
GHz ones for the a(f n i), but /„; can be significantly biased, 
mainly for the equilateral and the flat shape; on the other 
hand, at 225 GHz we observe an increase in the /„; error- 
bars. If we move far from this frequency range the capability 
of an accurate f n i measurement fast degrades, especially if 
only temperature is used. 

When polarization is considered, error-bars in the /„; 
estimate are reduced by at least a factor 2 for the ideal 
case (in agreemen t with estimates in literature; see e.g. 
iBabich et al.l 120041 ') and they are not affected by point 
sources in all the frequency range 75-225 GHz. Moreover, 
polarization can also help to significantly reduce the bias in- 
duced by unresolved sources. As example, for the equilateral 
shape, i.e. the shape most affected by point sources, the bias 
is a factor from 3 to 10 lower than from temperature only 
analysis. We can say therefore that polarization allows to 
expand the best range of frequencies where /„; can be safely 
determined. 

Discrepancies in the far -IR source bi s pectr um are ob- 
served when compared with lLacasa et all ([2012). They are 
attributed to the different ways to model the far-IR sources 
(see Section[S}. This dependency might affect the bias for the 
considered shapes especially at high frequencies (see Section 

E). 

Finally, as expected, the bias induced by point sources 
is strongly related to the flux density limit down to which 
sources are detected and removed/masked. Our results seem 
to point out that S c = 0.3 Jy is low enough to have unbi- 
ased /„; values at least at 165 and 225 GHz for a COrE-like 
mission. At different frequencies, lower flux density limits 
are required. On the other hand, /„; uncertainty is weakly 
affected by the choice of S c because it is dominated by the in- 
trinsic f n i variance (i.e. independent of contaminants), that 
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Figure 6. From left to right: the relative bias A(/ n ;)/o-(/ 7l ;) and the / n ; uncertainty cr(/ n ;) produced by unresolved point sources as a 
function of tmax for the ideal case. We consider the flux density limit S c =0.3 Jy and the following frequencies: 75, 105, 165, 225 GHz 
and the configuration without point sources (black line in right panels). From top to bottom, we plot results for the local, equilateral, 
orthogonal and flat f„i shapes. The solid lines take into account only temperature, while the long dashed lines include temperature and 
polarization. 
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Figure 7. As in Fig. [6] but for a COrE-like mission. 
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Table 4. The expected uncertainty <r(/ ni ), bias A/ n ; and relative bias A/ n ;/<r(/ n ;) at l m ax = 2048 from temperature only and 
temperature plus polarization for a COrE-like mission. S c = 0.3 Jy is adopted. 





case 


T 


T 


T 


T 


T 


T 


T+E 


T+E 


T+E 


T+E 


T+E 


T+E 




Freq. 


45 


75 


105 


165 


225 


375 


45 


75 


105 


165 


225 


375 




<*■(/«! ) 


9.7 


5.9 


4.7 


4.2 


5.3 


29.8 


6.0 


3.5 


2.7 


2.3 


2.5 


20.6 


local 




16.0 


4.3 


1.8 


0.4 


0.7 


25.1 


1.9 


-0.8 


-0.7 


-0.1 


-0.2 


8.2 




A/„i/cr(/ n! ) 


1.7 


0.7 


0.4 


0.1 


0.1 


0.8 


0.3 


-0.2 


-0.3 


-0.1 


-0.1 


0.1 




f(/nl) 


84.3 


65.2 


58.1 


56.0 


65.2 


232.6 


41.1 


27.5 


23.5 


21.0 


22.2 


147.2 


equilateral 


A/ ni 


421.9 


202.4 


102.7 


17.5 


49.0 


935.9 


66.2 


28.9 


13.5 


1.9 


4.6 


343.8 




Af nl /a(f nl ) 


5.0 


3.1 


1.8 


0.3 


0.8 


4.0 


1.6 


1.1 


0.6 


0.1 


0.2 


2.3 




f(/nl) 


49.7 


33.8 


28.4 


26.9 


32.0 


139.1 


21.7 


13.8 


11.4 


10.1 


10.6 


82.0 


orthogonal 


A/ ni 


39.9 


12.9 


6.5 


-0.1 


9.7 


367.8 


12.1 


23.8 


16.1 


2.8 


1.4 


130.2 




Af nl /a(f nl ) 


0.8 


0.4 


0.2 


-0.0 


0.3 


2.6 


0.6 


1.7 


1.4 


0.3 


0.4 


1.6 




v(fnl) 


93.8 


61.2 


50.9 


48.1 


57.7 


285.1 


39.1 


24.5 


20.2 


17.9 


18.9 


148.1 


flat 


A/ ni 


190.3 


68.1 


29.0 


6.5 


3.4 


-68.8 


10.4 


-26.3 


-20.4 


-3.7 


-5.2 


-38.4 




A/„,M/ n i) 


2.0 


1.1 


0.6 


0.1 


0.1 


-0.2 


0.3 


-1.1 


-1.0 


-0.2 


-0.3 


-0.3 



Table 5. The expected uncertainty <r(/ n ;), bias A/ n ; and relative bias A/ n ; /<r(f n i) at lmax = 2048 from temperature only and 
temperature plus polarization for COrE-like mission. S c = 1 Jy is adopted. 





case 


T 


T 


T 


T 


T 


T 


T+E 


T+E 


T+E 


T+E 


T+E 


T+E 




Freq. 


45 


75 


105 


165 


225 


375 


45 


75 


105 


165 


225 


375 




v(fnl) 


11.1 


6.4 


5.0 


4.4 


5.3 


29.9 


6.6 


3.7 


2.9 


2.4 


2.5 


20.6 


local 


A/„i 


130.7 


37.1 


15.6 


3.4 


2.3 


25.0 


15.3 


-7.0 


-6.0 


-1.2 


-0.6 


8.1 




A/„i/cr(/„i) 


11.8 


5.8 


3.1 


0.8 


0.4 


0.8 


2.3 


-1.9 


-2.1 


-0.5 


-0.2 


0.4 




v(fnl) 


92.4 


68.9 


60.6 


57.1 


65.7 


230.0 


43.7 


28.4 


24.0 


21.2 


22.3 


146.3 


equilateral 


A/ n! 


3428.4 


1731.9 


916.2 


175.0 


135.3 


933.6 


515.4 


235.0 


117.5 


19.5 


12.6 


346.7 




A/„j/(7(/„i) 


37.1 


25.1 


15.1 


3.1 


2.1 


4.1 


11.8 


8.3 


4.9 


0.9 


0.6 


2.4 




f(/ni) 


55.3 


36.2 


30.0 


27.5 


32.3 


138.7 


23.3 


14.3 


11.7 


10.2 


10.7 


81.8 


orthogonal 


A/ n ; 


553.1 


161.5 


79.9 


4.4 


16.4 


374.0 


110.3 


202.6 


136.9 


25.4 


12.2 


132.5 




A/ n j/<j(/ ni ) 


10.0 


4.5 


2.7 


0.2 


0.5 


2.7 


4.7 


14.2 


11.7 


2.5 


1.1 


1.6 




<r(fnl) 


105.8 


65.8 


53.9 


49.3 


58.3 


284.9 


42.0 


25.5 


20.8 


18.1 


19.0 


147.7 


flat 


A/ n ; 


1232.1 


523.1 


232.3 


58.3 


26.5 


-72.1 


59.4 


-227.2 


-172.4 


-33.0 


-14.7 


-39.3 




A/ ni /<r(/ ni ) 


11.6 


7.9 


4.3 


1.2 


0.5 


-0.3 


1.4 


-8.9 


-8.3 


-1.8 


-0.8 


-0.3 



decreases when increasing l mm . We want to stress the im- 
portance of high-resolution (i.e., FWHM<10 arcmin) and 
polarization-sensitive instruments, in order to have signifi- 
cant improvements in the precision of primordial /„; in fu- 
ture experiments beyond Planck. 
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